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High harmonic generation due to the interaction of a short ultra relativistic laser pulse with overdense 
plasma is studied analytically and numerically. On the basis of the ultra relativistic similarity theory 
we show that the high harmonic spectrum is universal, i.e. it does not depend on the interaction 
details. The spectrum includes the power law part /„ oc n -8 ' 3 for n < \/8a7n lax , followed by 
exponential decay. Here 7 max is the largest relativistic 7-factor of the plasma surface and a is 
the second derivative of the surface velocity at this moment. The high harmonic cutoff at oc 7^ ax 
is parametrically larger than the 47^ ax predicted by the "oscillating mirror" model based on the 
Doppler effect. The cornerstone of our theory is the new physical phenomenon: spikes in the 
relativistic 7-factor of the plasma surface. These spikes define the high harmonic spectrum and lead 
to attosecond pulses in the reflected radiation. 

PACS numbers: 42.65.Ky, 52.27.Ny, 52.38.Ph 
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I. INTRODUCTION 

High harmonic generation (HHG) from relativistically 
intense laser pulses interacting with solid targets has been 
identified as a promising way to generate bright ultra 
short bursts of X-rays [j], y, Q • For the first time this 
spectacular phenomenon was observed with nanosecond 
pulses of long wavelength (10.6 jum) CO2 laser light J4J. 

Short after the experimental observation in 1981, 
Bezzerides et al studied the problem of harmonic light 
emission theoretically [5| . Their approach based on non- 
relativistic equations of motion and hydrodynamic ap- 
proximation for the plasma predicted a cutoff of the har- 
monic spectrum at the plasma frequency. 

10 years later, in 1993, a new approach to the in- 
teraction of an ultra short, relativistically strong laser 
pulse with overdense plasma was proposed by Bulanov 
et al [6J. They "interpreted the harmonic generation as 
due to the Doppler effect produced by a reflecting charge 
sheet, formed in a narrow region at the plasma bound- 
ary, oscillating under the action of the laser pulse" [g- 
The " oscillating mirror" model predicts a cutoff harmonic 
number of 47^ ax , where 7 max is the maximal 7-factor of 
the mirror. 

At the beginning of 1996 numerical results of particle- 
in-cell simulations of the harmonic generation by fem- 
tosecond laser-solid interaction were presented by P. Gib- 
bon 0| ■ He demonstrated numerically that the high har- 
monic spectrum goes well beyond the cutoff predicted in 
|5| and also presented a numerical fit for the spectrum, 
which approximated the intensity of the n-th harmonic 
as /„ oc n~ 5 . At about the same time the laser-overdcnsc 
plasma interaction was also studied by Lichters et al |8j . 

The same year the analytical work by von der Lindc 
and Rzazewski Jj appeared. The authors used the "os- 
cillating mirror" model and approximated the oscillatory 
motion of the mirror as a sin-function of time without 
analysis of the applicability of this approximation. With 



the explicit form of the mirror motion an analytical for- 
mula for the harmonic spectrum was obtained. 

The first try to describe analytically the high harmon- 
ics generated at the boundary of overdense plasma by a 
short ultra intense laser pulse in a universal way that does 
not rely on an explicit formula for the exact plasma mir- 
ror motion was made in [lj . This work proposed for the 
first time the idea of universality of the harmonic spec- 
trum. In an attempt to advocate the "oscillating mirror" 
model and the role of the relativistic Doppler effect the 
authors of |l| made use of the steepest descent method 
and estimated the harmonic spectrum as /„ oc n~ 5 > 2 . 
However they did not discuss the area of applicability of 
the steepest descent method for the specific problem. We 
are showing in this article that applying the straightfor- 
ward steepest descent method is equivalent to applying 
the "oscillating mirror" model that proves to be insuffi- 
cient for the treatment of laser-overdense plasma inter- 
action. In the present work we revise [j| and clarify the 
physical picture of HHG. 

For the last couple of years the interest in the process 
of high harmonic generation from plasmas has enjoyed a 
revival thanks to the increasing interest in attoscience. 
The recent impressive progress in the physics of attosec- 
ond X-ray pulses [l(j triggers the fascinating question 
whether a range of even shorter ones is achievable with 
the contemporary or just coming experimental technol- 
ogy. 

L. Plaja et al 12| were the first to realize that the sim- 
ple concept of the " oscillating mirror" gives an opportu- 
nity to produce extremely short pulses and presented a 
numerical proof that the radiation generated by oscillat- 
ing plasma surfaces comes in the form of subfemtosecond 
pulses. For the first time the idea to use the plasma har- 
monics for the generation of subattosecond pulses (zep- 
tosecond range) was announced in [l|. 

The present article improves the analytical results of [lj 
and studies in detail the ultra relativistic plasma surface 
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FIG. 1: Geometry of the problem. The laser pulse is moving 
towards the overdense plasma slab, x is perpendicular to the 
surface, y and z are parallel to it. 



motion. We show that the high harmonics are generated 
due to sharp spikes in the relativistic 7— factor of the 
plasma surface. This new physical phenomenon leads to 
the spectrum cutoff at the harmonic number n cu toff 



"cutoff 



8a 7n 



(1) 



The cutoff QJ is much higher than the 47^ ax -cutoff pre- 
dicted by the simple "oscillating mirror" model. Here 
7max is the maximum 7— factor of the surface and a is 
a numerical factor of the order of unity, related to the 
plasma surface acceleration. 

Our analysis significantly exploits the relativistic 
plasma similarity theory |15|. which was developed after 
the work [lj had been published. The similarity theory 
enables us to rectify our previous results and to present 
them in a straightforward and clear way. 

In this article we first discuss the physical picture of 
high harmonic generation. The production of high har- 
monics is attributed to the new physical effect of the 
relativistic spikes. Then, we develop analytical theory 
describing the spectrum of the high harmonics and show 
that this spectrum is universal with slow power law decay. 
Finally, the theory is confirmed by direct particle-in-cell 
(PIC) simulation results. 



II. PHYSICAL PICTURE OF HHG AT 
OVERDENSE PLASMA BOUNDARY 

In this Section we state the problem of high harmonic 
generation at the boundary of overdense plasma and 
qualitatively describe its main features, which will find 
their analytical and numerical confirmation in what fol- 
lows. 

Let us consider a short laser pulse of ultra relativistic 
intensity, interacting with the sharp surface of an over- 
dense plasma slab (see Fig.[Q. 

We assume that the incident laser pulse is short, so 
that we can neglect the slow ion dynamics and consider 
the electron motion only. The electrons are driven by the 
laser light pressure, a restoring electrostatic force comes 
from the ions. As a consequence, the plasma surface os- 
cillates and the electrons gain a normal momentum com- 
ponent. 




FIG. 2: a) Electron momentum component parallel to the 
surface as a function of time b) The velocity of the plasma 
surface v s is a smooth function of time, unlike the 7-factor of 
the surface c) 



Since the plasma is overdense, the incident electromag- 
netic wave is not able to penetrate it. This means that 
there is an electric current along the plasma surface. For 
this reason, the momenta of electrons in the skin layer 
have, apart from the components normal to the plasma 
surface, also tangential components. 

According to the relativistic similarity theory |15J| . 
both the normal and tangential components are of the 
order of the dimcnsionless electromagnetic potential arj. 
Consequently, the actual electron momenta make a finite 
angle with the plasma surface for most of the times. 

Since we consider a laser pulse of ultra relativistic in- 
tensity, the motion of the electrons is ultra relativistic. In 
other words, their velocities are approximately c. Though 
the motion of the plasma surface is qualitatively differ- 
ent: its velocity v s is not ultra relativistic for most of the 
times but smoothly approaches c only when the tangen- 
tial electron momentum vanishes (see Fig. 12b). 

The 7-factor of the surface 7 S also shows specific be- 
havior. It has sharp peaks at those times for which 
the velocity of the surface approaches c (see Fig. |2t). 
Thus, while the velocity function v s is characterized by 
its smoothness, the distinctive features of 7 S are its quasi- 
singularities. 

When v s reaches its maximum and 7 S has a sharp peak, 
high harmonics of the incident wave are generated and 
can be seen in the reflected radiation. Physically this 



means that the high harmonics are due to the collective 
motion of bunches of fast electrons moving towards the 
laser pulse 13]. 

These harmonics have two very important properties. 
First, their spectrum is universal. The exact motion of 
the plasma surface can be very complicated, since it is 
affected by the shape of the laser pulse and can differ 
for different plasmas. Yet the qualitative behavior of v s 
and 7 S is universal, and since it governs the HHG, the 
spectrum of the high harmonics does not depend on the 
particular surface motion. 




FIG. 3: The universal high harmonic spectrum contains power 
law decay and exponential decay (plotted in log-log scale). 

We show below that the high harmonic spectrum con- 
tains two qualitatively different parts: power law decay 
and exponential decay (see Fig.[3J). In the power law part 
the spectrum decays as 



In OC 1/n 8 / 3 , 



(2) 



up to a critical harmonics number that scales as 7^ ax , 
where /„ is the intensity of the nth harmonic (see Sec- 
tion 0). Here 7 max is the maximal 7- factor of the point, 
where the component of the electric field tangential to 
the surface vanishes (see Section llVJ) . 

The second important feature of the high harmonics is 
that they are phase locked. This observation is of partic- 
ular value, since it allows for the_generation of attosecond 
and even subattosecond pulses 



III. ULTRA RELATIVISTIC SIMILARITY AND 
THE PLASMA SURFACE MOTION 



The analytical theory presented in this work is based 
on the similarity theory developed in |l5j for collisionless 
ultra relativistic laser-plasma regime and is valid both for 
under- and overdense plasmas. 

The ultra relativistic similarity theory states that when 
the dimensionless laser vector potential ao = eAo/mc 2 
is large (aft ^> 1) the plasma electron dynamics does not 



depend on ao and the plasma electron density N e sep- 
arately. Instead they merge in the single dimensionless 
similarity parameter S defined by 



S = 



N e 
a N c 



(3) 



where N c = ojQm/4ire 2 is the critical electron density for 
the incident laser pulse with amplitude ao and carrier 
frequency loq. 

In other words, when the plasma density N e and 
the laser amplitude ao change simultaneously, so that 
S = N e /aoN c = const, the laser-plasma dynamics re- 
mains similar. In particular, this basic ultra relativistic 
similarity means that for different interactions with the 
same S = const, the plasma electrons move along similar 
trajectories while their momenta p scale as 



p oc ao- 



(4) 



The S— similarity corresponds to a multiplicative 
transformation group of the Vlasov-Maxwell equations, 
which appears in the ultra relativistic regime. The sim- 
ilarity is valid for arbitrary values of S. Physically the 
S'-parameter separates relativistically overdense plasmas 
(S 3> 1) from underdense ones (S <§C 1). 

To apply the key result Q of the similarity theory to 
the plasma surface motion, we rewrite J3J for the electron 
momentum components that are perpendicular p n and 
tangential p T to the plasma surface: 



Pn oc Oo, 



p T oc a . 



(5) 



This result is significant. It shows that when we increase 
the dimensionless vector potential ao of the incident wave 
while keeping the plasma overdense, so that S = const, 
both p„ and p T grow as ao- In other words the velocities 
of the skin layer electrons 



P? 



mlc 2 + p 2 



P? 



= c(l - O(a 2 )) (6) 



are about the speed of light almost at all times. Yet the 
relativistic 7-factor of the plasma surface j s {t') and its 
velocity (3 a {t') behave in a quite different way. To realize 
this key fact let us consider the electrons at the very 
boundary of the plasma. The scalings JHJ state that the 
momenta of these electrons can be represented as 



p„(i') = a P n (S,LU t') 
p T (t') = aoP T (S,ui t'), 



(7) 



where P„ and P T are universal functions, which depend 
on the pulse shape and the ^-parameter rather than on 



a or N e separately. Consequently, for f3 s (t') and J s (t') 
one obtains 



&(*') 



Pn(t') 



^/m 2 e c 2 +p 2 n (t')+p 2 (t') 



P*(f)+P?(f) 



0(V 



(8) 



7.(0 



Pn(*0 



; = Wl + ^H + °(«0 )■ (9) 

VI - 0?(*') V P r(*0 



One sees from (0 that when <Iq gets large, the rela- 
tivistic 7-factor of the electrons becomes large too and 
their velocities approach the velocity of light. However, 
the dynamics of the plasma boundary is significantly dif- 
ferent. For large aos the plasma boundary motion does 
not enter the ultra relativistic regime and its relativis- 
tic 7-factor 7 s (t') is generally of the order of unity. Yet 
there is one exception: if at the moment if it happens 
that P T (S,t') = 0, i.e. 



Pr(S,C 



(10) 



we have 



One might expect that the "oscillating mirror" model 
could describe the laser-plasma interaction in our prob- 
lem. Therefore we want to explain in detail why it is not 
the case and then present the derivation of the correct 
boundary condition. 

The "oscillating mirror" model implies that the tan- 
gential components of the vector potential are zero at 
the mirror surface. As a consequence, if the ideal mir- 
ror moves with 7» 1 towards the laser pulse with the 
electric field E\ and the duration r , the reflected elec- 
tric field will be E ie ^ ex 7 2 -E; and the pulse duration 
will be -r re fl ex t / j 2 . The energy of the reflected pulse 
would be then j 2 times higher than that of the incident 
one. However, since the plasma surface is driven by the 
same laser pulse, this scaling is energetically prohibited 
and consequently the plasma cannot be described as an 
" ideal mirror" . 

To derive the correct boundary condition, let us con- 
sider the tangential vector potential components of a 
laser pulse normally incident onto a overdense plasma 
slab. These components satisfy the equation 



1 d 2 A T {t,x) d 2 A T {t,x) 4tt 



or 2 



dx 2 



j(t,x), (13) 



where A T (t, x = — oo) = and j is the tangential plasma 
current density. Eq. 1)13(1 yields 



7.s = 



'p£ 



V^M 



m±c* 



ex Oq. 



(11) 



So the relativistic 7-factor of the boundary jumps to 
7 S (t'g) ex ao and the duration of the relativistic 7-factor 
spike can be estimated as 



At ex l/(aoWo). 



(12) 



For the velocity of the plasma boundary one finds anal- 
ogously that it smoothly approaches the velocity of light 
as/3 s (i') = (1 — 0(a~ 2 )). Fig. [^represents schematically 
this behavior. 

As we will see later, the 7 S spikes cause the generation 
of high harmonics in the form of ultra short laser pulses. 



+00 



A T (t,x) = 2n / 3(t,x,t',x') dt'dx'. 



(14) 



Here J (t, x, t', x') 
have defined 0. 



j (t',x') (9- — 0+), where we 
= Q(t-t'-\x-x'\/c) and 6+ = 
8 (t — if + (x — x 1 ) /c), using the Heaviside step- function 
©(£). Due to this choice of J the vector potential A T (i, x) 
satisfies both Eq. (|13|) an d the boundary condition at 
X = — 00 since 3(t, x — —oo,t',x') — 0. The tangential 
electric field is E r = —(l/c)dtA T (t,x). If we denote the 
position of the electron fluid surface by X(t) we have 



E T (t,X(t)) 



2tt 



E 



j(t + «e/ C ,A(t) + 0<15) 



IV. BOUNDARY CONDITION: ENERGY 

CONSERVATION AND THE APPARENT 

REFLECTION POINT 

In this Section we introduce the boundary condition 
describing the laser-overdense plasma interaction appeal- 
ing to physical arguments, just as it was previously done 
in 0,[l3|- Mathematicall y rig orous analysis of the bound- 
ary condition is given in 13]. However, for the purposes 
of the present work it is sufficient to treat this problem 
on a more intuitive basis [l| . 



where f = x' - X(t). 

Now one has to estimate the parameters characterizing 
the skin layer, i.e. the characteristic time r s of skin layer 
evolution (in the co-moving reference frame) and the skin 
layer thickness S. Since the plasma is driven by the light 
pressure, one expects that r s ex 1/luq- The estimation of S 
is more subtle. From the ultra relativistic similarity the- 
ory follows that S ex (c/ujq)S a , where S ^> 1 for strongly 
overdense plasmas and A is an exponent that has to be 
found analytically. In this work we do not discuss the 
exact value of A, but notice that this quantity does not 
depend on neither S, nor oq. On the other hand, the 



denser the plasma, the less the value of 5. This condition 
demands that A < and we get c/uj ~^> 5 for S 3> 1. 

If the characteristic time r s of the skin layer evolu- 
tion is long (ct s 3> 5), then we can use the Taylor ex- 
pansion j (t ± £/c, x 1 = X(t) + £) m j(t, a;') ± e, where 
e = (£/c)dtj(t, x 1 ), and substitute this expression into 
(I15JI . The zero order terms cancel each other and we 
get E T (t,X(t)) ex J p (S/ct) <C Ei, where J p ex cEi is 
the maximum plasma surface current. Thus, as long as 
the skin-layer is thin and the plasma surface current is 
limited, we can use the Leontovich boundary condition 



M 



e„xE(f,I(i)) = 0. 



(16) 



V. HIGH HARMONIC UNIVERSAL 
SPECTRUM 

According to Eq. (|l(jfl . the electric field of the reflected 
wave at the plasma surface is 



E r (t',X{t')) = -E i (t',X(t')), 



(19) 



where E,(i',A(t')) = -(l/c)d t >Ai(t',X(t')) is the in- 
cident laser field, t' is the reflection time. The one- 
dimensional wave equation translates signals in vacuum 
without change. Thus the reflected wave field at the ob- 
server position x and time t is E r (i, x) — —~Ei(t',X(t')). 
Setting x — at the observer position we find that the 
Fourier spectrum of the electric field E r (£, x = 0) is 



This condition has a straightforward relation to energy 
conservation. Indeed, if we consider the Poynting vector 



-f oo 



47T 



-E x B. 



(17) 



we notice that the boundary condition ljTU|) represents 
balance between the incoming and reflected electromag- 
netic energy flux at the boundary X(t). 

The boundary condition (|16fl allows for another inter- 
pretation. An external observer sees that the electro- 
magnetic radiation gets reflected at the point x ARP (£), 
where the normal component of the Poynting vector 
S„ = cE r x B r /47r = 0, implied by E r (ir ARP ) = 0. 
We call the point x Anr (t) the apparent reflection point 
(ARP). 

The actual location of the ARP can be easily found 
from the electromagnetic field distribution in front of the 
plasma surface. The incident laser field in vacuum runs in 
the negative x— direction, E l (x,t) = ~E l (x + ct), while the 
reflected field is translated backwards: E r (a;, t) — E r (x — 
ct). The tangential components of these fields interfere 
destructively at the ARP position x ARr (t), so that the 
implicit equation for the apparent reflection point a; ARP (t) 
is: 



K(*a 



cf)+E;(x ARP -ct) = 0. 



(18) 



We want to emphasize that Eq. I|18|l contains the elec- 
tromagnetic fields in vacuum. That is why the reflection 
point x ARP is apparent. The real interaction within the 
plasma skin layer can be very complex. Yet, the external 
observer, who has information about the radiation in vac- 
uum only, sees that E r = at x ARP . The ARP is located 
within the skin layer at the electron fluid surface, which 
is much shorter than the laser wavelength for overdense 
plasmas, for which the similarity parameter is S 3> 1. In 
this sense, the ARP is attached to the oscillating plasma 
surface. 



E r (fi) 



eV2ir 



Re[ia((ci' +X(t'))/cr ) x 



x exp(— iujQt' — iujoX(t')/c)] exp(— i£lt) dt, (20) 



where 



t' -X(t')/c=t 



(21) 



is the retardation relation Q ■ 

The fine structure of the spectrum of E r (i) depends 
on the particular surface motion X(t), which is defined 
by the complex laser- plasma interaction at the plasma 
surface. Previous theoretical works on high order har- 
monic generation from plasma surfaces 0, H, EJ H2 tried 
to approximate the function X(t) in order to evaluate 
the harmonic spectrum. For the first time analytical de- 
scription of the high harmonic intensity spectrum and the 
concept of universality were presented in (ij . This work 
proclaims the idea that the most important features of 
the high harmonic spectrum do not depend on the de- 
tailed structure of X(t). However, the article 0] treats 
the universal spectrum without any relation to similarity 
theory for ultra relativistic plasmas and, as a result, re- 
lies on the saddle-point method without proper analysis 
whether or not the areas contributing to the spectrum 
overlap [j], [fj] • In our analysis we overcome these short- 
comings of 0| . 

To find the spectrum, we notice that the investigation 
of E r (0) (|20|) is equivalent to the investigation of the 
function 



/(n) = /+(n) + /_(«), 



(22) 



where 



+oo 



f± == ± I g(r +x(r ))exp(±i(r +x(t )) -wit) dr. 

(23) 




FIG. 4: Surface dynamics and path integration in 125H . a) 
Velocity x'(r') of the plasma surface; x' n — (n — l)/(n + 1) 
are the saddle points corresponding to dO/dr' = 0. b) The 
integration path can be shifted below the real axis everywhere 
except in the neighborhoods of r' g (dashed regions) . 



every smooth function resembles a parabola near its ex- 
terna, the expansion of x'(t') is a quadratic function of 
(t' — r'). Simple integration leads to the following ex- 
pression for x(t') 



a(r') 
x(t>) = x(r' g ) + v Q (r' g )(T' r' g ) - -^L(r> r' g )\ (27) 

The Taylor expansion given by Eq. I|27[) has three im- 
portant properties related to its dependence on 5* and 
a®: 1) for S = const and ao — > +oo one finds that 
vq — > c; 2) for ao — > +oo, a depends only on the parame- 
ter S; 3) the expansion (|27J) is a good approximation for 
\t' — T g \ <C (2n/ujo)fi(S), where the function /i does not 
depend on ao- These three properties are mathematical 
statements of the physical picture described in Section UTI 
combined with the similarity theory developed in Section 
IIIII In other words, the properties of the expansion 127|) 
just mentioned are direct consequences of the physical 
picture presented in Fig. [21 

Substitution of Eq. l(2Tjl into f± (n) yields 



Here r = Lotf, t' = u>ot\ n = £1/lo$, x(t') = (uio/c)X(t') 
and g is a slowly varying function (\dg(r')/dT'\ <C 1), 
which is trivially related to a as 



Z (t' + x{t')) = 



-im e c 



a((ct' + X(f))/en). (24) 



2eV2?r 
Making use of Eq. (|2~TJl we re-write Eq. l(2*3ll as 



+oo 



f± = ± 9 {t' + x(T'))exp{iT'{-n±l) + 



+ ix{ T ')(n ± 1)) (1 - x' (r')) dr'. (25) 

We wish to examine the integral l|25|l for very large n. For 
this purpose, we notice that the derivative of the phase 



/±(n) = £/ ± (r», 



where the sum is over all times r' 



F{r'n) = ^ Ai 



(28) 



f+i T ' g i n ) = 9 {-r'g + x(Tgj) exp{iO + (T g ,n)) x 

xF(T' g7 n) (29) 

f-{ T ' g , n ) = -9(r g + x(r g ))exp(i@_(r' g ,n))x 

xF( Tg ,-n) (30) 

2 n-n cr (T' g ) \ 
a(T>yni^"\ri cr (T>)(a(Ti)ny/i) y ' 



6± = ±(T g +x(r' g ))+n(x(r' g )-T' g ), (32) 

and n cr — 2/(1 — vo). In (|31|l Ai is the well-known Airy- 
function, defined as 



e(r') = T '(-n ± 1) + x(T r ){n ± 1) 



(26) 



is negative everywhere except in the vicinity of t'„ = loqI' 
for which x'{T g ) « 1 (see Fig. QJi). 

The physical meaning of t' and the behavior of x(t') 
in the vicinity of these times is explained by Eq. (|10fl . 
Since the time derivative of O(r') is negative for all rs 
that are not too close to one of the r' we can shift the 
path over which we integrate to the lower half of the 
complex plane everywhere except in the neighborhoods 
of t' (see Fig.0J}). The contributions of the parts remote 
from the real axis are exponentially small. We can shift 
the path to the complex plane till the derivative equals 
zero or we find a singularity of the phase 0. 

To calculate the contributions of r' g s neighborhoods we 
can expand x'(t') near each of its maxima at r' Since 



Ai(x) = —j= \ c 
yvr J 



cos ( ux + — u J ) du. (33) 



Note that if x(t' +tt) = x(t') and g{r') = g(r' + tt), then 
/±(2n) = 0. 

Using equations l|28|) - (|32|l we can show analytically 
that the spectrum of radiation generated by the plasma 
is described by a universal formula. 

For the intensity of the nth harmonic we obtain 



/„ oc |exp(z6 + (n)) F (n) — exp (i0_(n)) F (— n)|" 



(34) 



where 



F(n) 



4^ 



Ai 



/on) 4 /3 
± (n) = ±00-7101, 



2 n 



(an)V3 



(35) 
(36) 



with the Airy function Ai(x) defined in 1331) and the crit- 
ical harmonic number n cr satisfying n cr = 47^ ax , where 
7max is the largest relativistic factor of the plasma bound- 
ary. 

Eq. I|34|) gives an exact formula for the high harmonic 
spectrum, which includes both power law and exponen- 
tial decay parts. Now we want to use different asymptotic 
representations of the Airy function in order to demon- 
strate explicitly these two quite different laws of high 
harmonic intensity decay. 



For n < 



n/n cr \ < (an) 1 / 3 ) 



wo 



can substitute the value of the Airy function at x = 
(Ai(0) = V^/(3 2/3 r(2/3)) = 0.629, Ai'(0) = 
-3 1 / 6 r(2/3)/(2 v ^F) = -0.459) in Eq. g3), and obtain 



oc 



,8/3 



sin Oq 



Ai'(0) 
Ai(0) 



B(n,Q ) 



(37) 



where 



B(n,@ ) 



2 sin 6o 



2/3 



2i cos 6o 



(ancr) 1 / 3 \n cr J (an) 1 / 3 

This means that the universal spectrum 

/„ oc 1/n 8 / 3 



(38) 



(39) 



is observed everywhere except for sinGo ~ 0, when the 
dominant term in the expansion is zero. For this partic- 
ular case, a higher order correction is important for 



(an) 1/3 tan0 o < 2 



Ai'(0) 



Ai(0) 



(40) 



and in this restricted frequency range the spectrum 



J„ oc 1 /n 10 / 3 



(41) 



should be used. 

At this point we want to explain the meaning of 8/3- 
spcctrum universality. Notice that since Eq. (|37|l de- 
pends on the phase 0, for moderate values of n cr the 
best power law fit of the high harmonic spectrum can be 
delivered by 



oc n 



(42) 



where 8/3 < p < 10/3. When n cr becomes really 
large, the majority of the harmonics does not satisfy 



the inequality l|40(l and the spectrum inevitably becomes 
I n oc n~ 8 / 3 . In other words, one should think of the 8/3- 
spectrum as the high intensity limit of the high harmonic 
spectrum. To find an analytical criterion for the 8/3- 
spectrum generation one can state that the condition l|40() 

has to be violated for the harmonics with n oc 
This means that 



In 



>78/3 



1 


Ai'(0) 
Ai(0) 


\/2tan0o 



(43) 



The formal expansion of the Airy function for n <C n cr , 
(an)}/ 3 < 1 leads to the spectrum I n oc 1/n 5 / 2 discussed 
in [l|. Meticulous analysis demonstrates that this case 
is irrelevant. Yet the version of the saddle point method 
used in [l[ ascribes the spectrum 1/n 5 / 2 to the whole area 
below the cutoff. This is because the areas around the 
saddle points contributing to the spectrum overlap signif- 
icantly and the approximation used in [l| is not accurate 
enough in this part of the spectrum. However, the differ- 
ence between the powers 8/3 and 5/2 is only 1/6 and it is 
hardly distinguishable numerically and experimentally. 

For n > ^faj&nll 2 (2|1 - n/n cr \ > (an) 1 / 3 ) Eq. (JHH) 
can be rewritten as 



1/2 



■ exp 



lev 7 ^, 



3a 1 / 2 3 / 2 



(44) 



It is interesting to notice that the approximation used 
in [jj also gives Eq. l|44|l for this area, i.e. here the 
overlapping of the contributing areas is really negligible. 



VI. ULTRA SHORT PULSE DURATION 

For ultra short pulse generation not only the ampli- 
tudes but also the harmonic phases are of importance. 
The calculations presented above show that all harmonics 
are phase locked. This means that after proper filtering 
they can produce a pulse of duration T, such that 



Toe 



1 



1 



oc 



UJQ 



ll 2 Vo-% 



3 ' 
max 



(45) 



where w is the frequency of the fundamental wave. 

Eq. (|45|) presents a new result. Notice that the plasma 
boundary never moves with a relativistic 7-factor larger 
than 7 m ax- Consequently, the frequency of a photon re- 
flected from this boundary due to the relativistic Doppler 
effect does not exceed 47 max cjo = n cr ujQ. How can a pulse 
with duration T given by Eq. (|45|l be produced then? 

Mathematically this can be understood looking at the 
properties of the Airy function. Ai(x) changes its behav- 
ior from oscillatory for x < with |ar| 3> 1 



Ai(x) 



1 



1 1/4 



?w 3 / 2 + - 



(46) 
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to exponentially decaying for i>1 



Ai(x) 



2x x l i 



exp 



r 3/2 



(47) 



The point x — corresponds to n — n cr . However, the 
exponential is so small for n < yo/8n C r that the power 
law decay dominates over the exponent. As a result, the 
power law spectrum goes well beyond the threshold u)on cr 
predicted by oversimplified models considering reflection 
from a moving surface. 

The reasoning just presented explains the mathemat- 
ical origin of the ^JaJ&n C r -cutoff. In the next section 
we give a simple physical interpretation of this result and 
reveal its relation to the relativistic 7-factor spikes. We 
also take a closer look at the mathematical nature of the 
cutoff and its relation to the failure of the standard saddle 
point approach due to saddle point overlapping. 



VII. RELATIVISTIC SPIKES AND CUT-OFF 
OF THE HIGH HARMONIC SPECTRUM 

The 7 3 — scaling of the spectrum cutoff 0J is readily 
understood using the relativistic 7-factor spikes of the 
plasma surface motion. Indeed, the plasma surface ve- 
locity in the vicinity of a maximum can be approximated 

as 



V{t) = V {tg) - Otwl{t - tgf 



(48) 



(see Section 0. Consequently, for the surface 7-factor 
during a relativistic spike we find 



lit) 



V 1 +lm^a^(t-t g ) 27 



(49) 



where 7 max = 1/yl - -u 2 (t s )/c 2 . Eq. l{4*9"|) shows that 
the highest harmonics are generated over the time inter- 
val 



what physically explains the origin of the high frequency 
cutoff jl) . This cutoff should be compared with the one 
predicted by the "oscillating mirror" model: 47^ ax . It 
differs parametrically from the correct cutoff tx v8a7^ lax , 
which is due to the relativistic 7-factor spikes. 

As it has been shown, the "oscillating mirror" model 
gives the incorrect formula for the spectrum cutoff be- 
cause it does not include the relativistic 7-factor spikes. 
Mathematically, this failure " oscillating mirror" model is 
related to the saddle point overlapping. Eq. i|2ti|) defines 
the saddle points 



Mr') _ 1± 1 
dr' n 



(53) 



the vicinity of which determines the value of the integral 
describing the amplitude of the n-th harmonic, n ^$> 1, 
see Fig.Q] Eq. (J57J) yields that the set of Eqs. J53J has real 
(not imaginary) solutions only for n < n cr = 47^ ax . For 
larger n all of the saddle points of (|53[1 have an imaginary 
part. 

Let us now apply the standard saddle point method 
to the problem without investigating whether the condi- 
tions for its applicability are met (it is clear that these 
conditions are violated at least for n ~ n cr , for which 
two saddle points coincide) . This approximation predicts 
that for 1 -c n <C n cr the spectrum decays as 1/n 5 ' 2 . For 
n > n cr this approach restores Eq. (|47|l . 

As it was mentioned above the spectrum 1/n 5 / 2 occurs 
formally in the limit (an) 1 ' 3 < 1. In other words, if a 
were small, this spectrum would be observed. However, 
since a depends on S rather than ag, this spectrum cor- 
responds to harmonics with small numbers only and is 
hardly of any practical interest. On the other hand, one 
can notice that a describes the plasma surface accelera- 
tion at the maximum of the velocity. This means that the 
limit of small a (and the spectrum 1/n 5 / 2 ) describes the 
limit of the surface moving without acceleration. How- 
ever, this limit is of little interest for large ag. 
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wo \/ajn 



For the whole time interval At the relativistic spike moves 
with ultra relativistic velocity in the direction of the emit- 
ted radiation. For this reason, the spatial length L of the 
high harmonic pulse produced is 



L oc (c — vo(t g ))At oc 



1 



^0 V^Tmax 

A pulse of such duration contains frequencies 



n oc - oc woV^nd^, 



VIII. CUTOFF AND THE STRUCTURE OF 

(50) FILTERED PULSES 

As we have seen, the relativistic plasma harmonics are 
phase locked and can be used to generate ultra short 
pulses. However, to extract these ultra short pulses, one 
has to remove the lower harmonics. The high energy 
cutoff Q of the power law spectrum defines the shortest 
pulse duration that can be achieved this way. 

(51) Let us apply a high-frequency filter that suppresses all 
harmonics with frequencies below £lf and study how the 
relative position of the fi/ and the spectrum cutoff affects 
the duration of the resulting (sub)attosecond pulses. 

According to Eq. Q34JI. the electric field of the pulse 
(5^) after the filtration is 



IX. SPECTRUM MODULATIONS 



+ oo 



E ex Re I (exp (i6+(n)) F{n) — exp(i0_(n)) x 

Slf/uJo 

x.F(— n)) exp(ini) (in (54) 

The structure of the filtered pulses depends on where 
we set the filter threshold 17/. In the case 1 <§; 17//cjo <C 
^a/8n c ' r , we use Eq. (|55|) and rewrite Eq. IJ51|) as 
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EocRe 
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TITS - 
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n^/w 
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1/3 



-^- Reexp{ittft-i®i)P(ttft), (55) 



where the function P 



P{x) 



+00 



exp(iyx) 



y 



4/3 



dy 



(56) 



gives the slow envelope of the pulse. 

It follows from the expression 1|55J) that the electric 
field of the filtered pulse decreases very slowly with the 

— 1/3 

filter threshold as u. ' . The pulse duration decreases 
as 1/17/ . At the same time, the fundamental frequency 
of the pulse is 17/. Therefore the pulse is hollow when 



Qf/ujo <S Va874axi i- e - its envelope is not filled with 
electric field oscillations. One possible application of 
these pulses is to study atom excitation by means of a 
single strong kick. 

The pulses structure changes differ significantly when 
the filter threshold is placed above the spectrum cutoff. 

For Qf/uo > Va87max we use E q s - <EH> and CD to 
obtain 



E ex 



up 



3/2 



exp 



xRe- 



8V2CL f 
exp(iftft — i@i) 



8a/2/ (3^ 



3/2 



iujot 



(57) 



Th amplitude of these pulses decreases fast when 17/ 
grows. However, the pulse duration ex l/v^7max does 
not depend on 17/. Since the fundamental frequency 
of the pulse grows as 17/, the pulses obtained with an 
above-cutoff filter are filled with electric field oscillations. 
Therefore these pulses are suitable to study the resonance 
excitation of ion and atom levels. 

The minimal duration of the pulse obtained by cutting- 
off low order harmonics is defined by the spectrum cut- 
off ex \/a87m ax - Physically, this result is the conse- 
quence of the ultra relativistic spikes in the plasma sur- 
face 7— factor. 



In this section we discuss the harmonic phases and 
show how the interference of harmonics produced by dif- 
ferent 7— spikes can lead to spectrum modulations. 

Eq. (|28[) obtained in Section [V] allows for straight- 
forward physical interpretation. The sum /_(r',n) + 
f + (r g ,n) gives the contribution of the r^-spike of the 
surface relativistic factor to the harmonic spectrum (see 
Fig-EJ- Therefore the phase of the nth harmonic 4> n due 
to the T' g s spike is given by 



tan(0 n (rp + n@i(T' g )) tan (6 ) = 

= l-F{r' g ,-n)/F{T' g ,n) 
l + F(T'-n)/F(r'n)- 



(58) 



where Qq(t') = r' + x(t') is the phase of the incident 
laser pulse at the time r g and @i(t') — T g —x(T g ) = r(T' g ), 
the time of the observation, i.e. the symbol Oi is redun- 
dant, yet we keep it for the sake of notation coherency, 
-n) « F(T' p ,n) for n > 1, Eq. JSSJ gives 



Since F{t' 
J- 



rise to ^> n {T g ) ~ — n@i(r g ). This means that each spike 
radiates phase locked high harmonics that can be used 
for ultra short pulse production. Another consequence of 
F(t'~ti) w F(t'ti) is that the amplitude of the har- 
monics is proportional to sin {@o(t')). 

The mechanism presented in Fig. [21 has another very 
interesting consequence. Each harmonic is generated 
due to several spikes. These spikes contribute to l|28|) 
with different phase multipliers. This phase difference 
leads to modulations in the high harmonic spectrum. 
As an example we consider the interference between the 
harmonics produced by two different spikes in detail. 
The phase shift between the contributions from differ- 
ent spikes is <f> n (T gi ) - Kij'g^) = @o(r gi ) - ©o(^ 2 ) - 
n(Qi(r gi ) — Qi(T g2 )) and can be large for large n. Since 

for n <C min 



^a(r gi )/8nU 2 (r' gi ), ^a(r g2 )/8n 3 c l 2 (r' g2 



only Ai(0) enters /±(t' ), the values of the contribu- 
tions from the r' spikes for this harmonic range do not 
depend on n cr (t' ) . As a result the modulations in this 
harmonic range depend only on the parameter S. 

To recapitulate, the non-trivial plasma motion produc- 
ing more than one 7-spike per oscillation period is the 
cause of the spectrum modulation. 



X. NUMERICAL RESULTS 

In order to check our analytical results, we have 
performed a number of Id PIC simulations with the 
Id particle-in-cell code VLPL 18]. In all simula- 
tions a laser pulse with the Gaussian envelope a = 
ao exp (— i 2 /r|) cos (woi), duration cj tl = 4-7T and di- 
mensionless vector potential a = 20 was incident onto a 
plasma layer with a step density profile. 
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A. Apparent Reflection Point 

First, we study the oscillatory motion of the plasma 
and the dynamics of the apparent reflection point de- 
fined by the boundary condition (|l(j|l [20| . The plasma 
slab is initially positioned between xr = — 1.5A and 
Xl = — 3.9A, where A = 2tt/ujq is the laser wavelength. 
The laser pulse has the amplitude ao = 20. The plasma 
density is N e /N c = 90 {S = 4.5). 
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FIG. 5: ID PIC simulation results for the parameters ao = 20 
and N e = 90-/V c . a) Oscillatory motion of the point x A Rp(t) 
where E T (x(t)) = O.b) Velocity i> A Rp(t) = dx AIl p(t)/dt; only 
the positive velocities are shown, since they correspond to 
motion towards the laser pulse in the geometry of this simu- 
lation. Notice that the ARP velocity is a smooth function, c) 
The corresponding 7— factor 7arp(£) = 1/yl — v A Rp(t) 2 /c 2 
contains sharp spikes, which coincide with the velocity ex- 
trema. 

At every time step, the incident and the reflected fields 
are recorded at x — (the position of the " external ob- 
server"). Being solutions of the wave equation in vac- 
uum, these fields can be easily chased to arbitrary x and 
t. To find the ARP position a? ABP , we solve numerically 
equation (|18fl . The trajectory of x ARP (t) obtained in this 
simulation is presented in Fig.[SJi. One can clearly see the 



oscillatory motion of the point x AKP (t). The equilibrium 
position is displaced from the initial plasma boundary 
position xr due to the mean laser light pressure. 

Since only the ARP motion towards the laser pulse is 
of importance for the high harmonic generation, we cut 
out the negative ARP velocities v A¥iP {t) = dx ABP (t)/dt 
and calculate only the positive ones (Fig. |5p) . The corre- 
sponding 7- factor 7arp(£) = 1/yl — v A Rp(t) 2 /c 2 is pre- 
sented in Fig. 0:. Notice that the ARP velocity is a 
smooth function. At the same time, the 7-factor 7 A rp(£) 
contains sharp spikes, which coincide with the velocity 
extrema. These numerical results confirm the predictions 
of the ultra relativistic similarity theory, which were pre- 
sented in Section ITTT1 



B. High harmonic spectrum 

For the same laser-plasma parameters (ao = 20, N e = 
90N C ) the spectrum of high harmonic radiation is pre- 
sented in Fig. El The maximum 7— factor of the ap- 
parent reflection point in this numerical simulation is 
7max ~ 3.3 (compare with Fig. |5J). Consequently the 
maximal harmonic number predicted by the "oscillating 
mirror" model lies at ^ 2 lax ~ 40, while the harmonic 
cutoff predicted by the relativistic spikes is about 100. 
Fig.[£j]clearly demonstrates that there is no change of the 
spectrum behavior at 47^ a£C , while steeper decay takes 
place above 100, as predicted by our theory. Also, the 
spectral intensity modulations discussed in SectioulVland 



spectri 



are observed. 




co/co„ 



FIG. 6: Spectrum of high harmonics obtained numerically for 
the case of ao = 20 and iV e = 90A^ C , corresponding to S — 4.5 
and jmax ~ 3.3. Assuming a « 1, the cutoff (1) is expected 
at n ~ 100. This analytically predicted cutoff is marked by 
the dashed line. 

To be able to make a real statement about the power 
in the power law decay of the spectrum we need more 
harmonics in order to satisfy the condition of universal 
8/3-spectrum formation l|43ll . For this reason we made 
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FIG. 7: Spectra of the reflected radiation for the laser am- 
plitude do = 20 and the plasma density N e = 307V cr . The 
broken line marks the universal scaling I oc cj~ a ' 3 . 



the simulation with parameters ao = 20 and N e — 30N C , 
which roughly corresponds to solid hydrogen or liquid 
helium. The reflected radiation spectrum obtained for 
these parameters is shown in Fig. 0in log- log scale. The 
power law spectrum /„ oc 1/n 8 ' 3 is clearly seen here, thus 
confirming the analytical results of Section Ivl 



C. Subattosecond pulses 

Let us take a closer look at Fig. \7\ The power law 
spectrum extends at least till the harmonic number 2000, 
and zeptosecond (lzs = 10~ 21 s) pulses can be generated. 
The temporal profile of the reflected radiation is shown 
in Fig. |S1 When no spectral filter is applied, Fig. [SJi, a 
train of attosecond pulses is observed [l^] • 
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FIG. 8: Zeptosecond pulse train: a) temporal structure of 
the reflected radiation; b) zeptosecond pulse train seen after 
spectral filtering; c) one of the zeptosecond pulses zoomed, its 
FWHM duration is about 300 zs. 



D. Filter threshold and the attosecond pulse 
structure 



The dependence of the short pulses on the position 
of the filter also can be studied numerically. We apply 
a filter with the filter function f(u) = 1 + tanh((w — 
ftf)/Au)). It passes through frequencies above f2/ and 
suppresses lower frequencies. We choose the similation 
case of laser vector potential ao = 20 and plasma density 
./V e = 90 N c . The spectrum of high harmonics is given in 
Fig. El 
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FIG. 9: Dependence of the pulses filling on the position of 
the sharp filter boundary for ao = 20 and N e = 90iV c and 
filter positions: a) fi/ = 20loq, Alo — 2ujo; b) Clf — 40wo, 

Auj = 2lu ; c) Q f = IOOcjo, Alo = 2uo Q ; d) tt f = 200lu , 
Alo = 2loo 

We zoom in to one of the pulses in the pulse train ob- 
tained and study how the shape of this one pulse changes 
with flf. Fig.|^lrepresents the pulse behavior for four dif- 
ferent positions of ilf . We measure the degree of fillness 
of the pulse by the number of field oscillations within the 
FWHM. One clearly sees that for filter threshold below 
the cutoff frequency, Fig. El a), b), the pulse is hollow. 
Notice that the case of Fig. ^jp corresponds to the cutoff 
frequency predicted by the "oscillating mirror" model. 
Only for filter threshold positions above the spectrum 
cutoff given by (JIJ the pulse becomes filled, Fig. |3;,d. 
These results confirm once again the real position of the 
harmonic cutoff. 



However, when we apply a spectral filter selecting har- 
monics above n — 300, a train of much shorter pulses 
is obtained, Fig. |Sp. Fig. IHJ; zooms in to one of these 
pulses. Its full width at half maximum is about 300 zs. 
At the same time its intensity normalized to the laser 
frequency is huge {eE zs /mcLo) 2 m 14. This corresponds 
to the intensity I zs w 2 x 10 19 W/cm 2 . 



XI. DISCUSSIONS 

In this work we have shown analytically and numeri- 
cally that the relativistic 7-factor spikes are the physical 
cause for high harmonic generation at the boundary of 
overdense plasma. It is important that the properties of 
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these spikes are universal and follow from the ultra rel- 
ativistic similarity theory. The universal physics of the 
relativistic 7-spikes inheres in the universality of the high 
harmonic spectrum. 

The spectrum of the high harmonics contains the 
power law part /„ ex I/71 8 / 3 , which goes till the cutoff at 
V / 8a7maa:- Here 7 maa; is the maximal 7-factor and a de- 
scribes the acceleration of the plasma boundary. This re- 
sult demonstrates that a naive " oscillating mirror" model 
is insufficient for correct treatment of high harmonic gen- 
eration at plasma boundaries. 

It is interesting to note though that if the plasma 



boundary moves without acceleration (a — > 0) our ap- 
proach restores the cutoff 4 7^^ following from the " os- 
cillating mirror" model and leads to I n ex 1/n 5 / 2 , yet this 
limit is irrelevant to laser-relativistic plasma interaction. 
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